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Abstract. The effect on particle physics processes by intense electromagnetic fields in the charge bunch collisions at 
future lepton colliders is considered. Since the charge bunch fields are tied to massive sources (the e + e“ charges), a 
reference frame is chosen in which the fields appear to be co-propagating. Solutions of the Dirac equation minimally 
coupled to the electromagnetic fields reasonably associated with two intense overlapping charge bunches are obtained 
and found to be a Volkov solution with respect to a null 4-vector whose 3-vector part lies in the common propagation 
direction. These solutions are used within the Furry interaction picture to calculate the beamstrahlung transition rate 
for electron radiation due to interaction with the electromagnetic fields of two colliding charge bunches. New analytic 
expressions are obtained and compared numerically with the beamstrahlung in the electromagnetic field of one charge 
bunch. The techniques developed will be applied to other collider physics processes in due course. 


1 Introduction 

Physics experiments in particle colliders involve the overlap of 
dense bunches of charged particles which interact to produce 
final states of interest. The treatment of the individual inter¬ 
actions which produce new particles is of course well known 
and thoroughly carried out in physics analysis. By comparison, 
the treatment of collective effects from the strong electromag¬ 
netic fields has been limited to beamstrahlung and background 
pair production processes only. Nevertheless, the overwhelm¬ 
ing bulk of physics processes take place in strong overlapping 
charge bunch fields. There is a lack of a general analysis in this 
respect and it is the goal of this paper to lay the groundwork of 
just such a general analysis. 

At linear colliders, beamstrahlung and coherent pair pro¬ 
duction have been considered to take place in a constant crossed 
electromagnetic field. Theoretically, the quasi-classical approx¬ 
imation (QCA) has been used to determine the transition rate 
mm. A method which does not rely on kinematic approx¬ 
imations and which takes into account the external field non- 
perturbatively, the Furry picture |4), was used to calculate the 
same process transition rates 0. The two methods produce 
identical results after the kinematic approximations of the QCA 
are accounted for 0. 

The equivalent photon approximation (EPA) has often been 
utilised to calculate the effect of charged particle fields on col¬ 
lider processes. The main processes producing e + e _ pairs are 
the Landau-Lifshitz and Bethe-Heitler processes. At linear col¬ 
liders these processes, along with the Breit-Wheeler process 
make up the incoherent pair production 0 . The applicability of 
the EPA depends on the size of the impact parameter compared 
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to the Compton wavelength and is not suitable for multipho¬ 
ton radiation 0- The inapplicability of the EPA also includes 
beamstrahlung, neutrino and heavy lepton pair production 0. 

In this paper, the effect of charge bunch fields on general 
processes of any order, without kinematic approximations, will 
be considered. The natural theoretical framework will be the 
Furry picture which requires solutions of the field equations 
embedded in the electromagnetic fields of the charge bunches. 
There appears to be, up to the present, no consideration of 
the case of overlapping, counter-propagating charged bunches. 
This is the particular problem to be addressed in this paper. 

Exact solutions of the fermion equations of motion with 
external electromagnetic field are known for plane wave fields, 
centrally symmetric fields and a few cases of a combination of 
fields ID- Exact solutions for spin 0 and spin 1 charged parti¬ 
cles embedded in plane wave fields are also known flOl . 

For the case of two plane waves, exact fermion solutions 
have previously been found for the restricted case of collinear 
fields with orthogonal vector potentials 111 II 121 , Volkov solu¬ 
tions of combinations of laser fields of different modes have 
been used in the study of photon emission and pair photo¬ 
production m. The case of two anti-collinear, circularly po¬ 
larised fields was considered by oa for the Klein-Gordon equa¬ 
tion, and by D3 for the Dirac equation. 

The case of overlapping charge bunch fields dealt with in 
this paper, will be considered in a reference frame in which 
the required computations are simplest. The Lorentz transfor¬ 
mation of photons reveals that two, free electromagnetic fields 
either co-propagate in all reference frames or they can be trans¬ 
formed to a frame in which they anti-propagate. The fields 
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associated with charge bunches, however, are tied to massive 
sources. It is always possible therefore to adopt a reference 
frame in which the charge bunches and their associated fields 
appear to co-propagate. 

In sectioned will be argued that it is not possible to assume 
that the vector potentials of the two charge bunch fields are 
orthogonal. Lepton colliders bring together bunches of oppo¬ 
site charge whose 3-potentials are generally anti-parallel. How¬ 
ever, the bunch collisions are not perfectly head-on, and the 
bunches distort during the collision. The 3-potentials of the two 
bunches, though co-planar, must otherwise be allowed to have 
general orientation in analytic calculations. 


The solutions sought in section [3] will be those of charged 
spin 1/2 particles in two co-propagating plane wave fields of 
otherwise general form. In such a combination of fields, a Volkov 
solution with respect to a null 4-vector whose unit 3-vector part 
lies in the common propagation direction, will be employed. In 
order to be certain this procedure is valid, an explicit proof is 


given in Appendix A 


In order to obtain commutation relations and quantisation 
of wavefunctions for use in the calculation of transition rates, 
orthogonality and completeness of the exact fermion solutions 
in two fields is necessary. The proof can be considered to be the 
same as that for the Volkov solution for a fermion in a single 
plane wave field. It can then be assumed that the commutation 
relations and quantisation within the Furry picture are canoni¬ 
cal. 


The explicit fermion solutions in two charge bunch fields 
will be applied directly to the calculation of the beamstrahlung 
transition rate in section [4] A new analytic form for this tran¬ 
sition rate will be obtained and compared to the usual beam¬ 
strahlung expression. 

This paper intends to lay the groundwork for a consider¬ 
ation of general collider phenomenology in the overlapping 
electromagnetic fields of colliding charge bunches. Given the 
strength of the fields expected at future linear colliders, this an¬ 
alytic work is considered essential, particularly for precision 
calculations. The next steps in this analytic program, following 
on from this paper, will be discussed in section [5] 

This paper works in natural units and the electromagnetic 
gauge is fixed by a combination of the Lorenz gauge and setting 
the time component of the 4-potential to zero, 

c = h = 47reo = 1, e = \]a 

0, A e = {0,A e ) (1) 


2 The external field of two colliding charge 
bunches 

Lepton colliders, from the point of view of the laboratory ref¬ 
erence frame, bring together two ultra relativistic, Gaussian 
charge bunches with a flat transverse profile. Since it is the aim 
of this paper to consider particle processes in the overlapping 
fields of both charge bunches, it is sensible to examine the col¬ 
lision in a reference frame that simplifies the calculation. 

By moving the point of view relativistic ally, it is possible to 
adopt reference frames in which either of the bunches appear 
to be at rest (the I\' and K" frames). By moving the observer 
further the K frame can be adopted in which both bunches ap¬ 
pear to be moving relativistically and collinearly (figure [lj. 


centre of mass frame 



K’ frame 

P’ 2 




K” frame 



Fig. 1. Apparent charge bunch motion for an observer in center of 
mass, K, K’ and K” reference frames. 



Fig. 2. Volume elements and associated vectors in the K frame. 


Superscripts, i/j ± , u ± will refer to the positive and negative 
energy solutions of the Dirac equation. Subscripts consisting 
of Greek letters will refer to 4-vector components in Cartesian 
coordinates and ± subscripts refer to linear combinations of 4- 
vectors and functions. Two or three component vectors will be 
represented by, x. 


In order to study particle physics processes in the fields of 
the bunches it is desirable to know the nature of the fields seen 
at an observer point O in the in the K frame. Infinitesimally 
small volume elements,Vi and V 2 , containing charge q\ and <72 
move with relativistic beta /3 i,/ 32 collinearly along the z-axis 
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(figure|2|. The vectors r\. r'j, in the K' and K" frames respec¬ 
tively, transform to transverse (ru, r 2 j_) and longitudinal (z) 
vectors in the K frame as, 

r’ 1 =r l _ L +ii{z-v 1 t)z, r 2 =r 2± +'y 2 (z -v 2 t)z ( 2 ) 

so that the electric fields in each element V\ and V 2 , within 
the K' and K" frames respectively, are 

E ! _ gir'i _ giri _L +yi(z-vit)z 

l r il 3 [|t-ix | 2 +7i(z - vit ) 2 ] S|/2 

■pH _ <n r 2 _ 92T-2X + 72(2 - V2t)z 

\ r 'i\ S [k 2 ±| 2 +7|(2 - v 2 t) 2 ] 3/2 

The transformation of the fields themselves to the K frame 
is achieved with the relations ( lfl 6 ) . section 11 . 10 ) 

E 1 = 7i £1 , 1 - i ^-/3i(/3i -E[)- B 1 =f3 1 xE 1 

E 2 = l2 E'i--^p 2 (p 2 -E'Z)- B 2 = (3 2 xE 2 (4) 

The combined electric and magnetic fields of the two vol¬ 
ume elements seen by the observer in the K frame is, 

E = E-\ + E'2 

_ 71 (gi r i± + (z - vit)z) 'y 2 (q 2 r 2 ± + (z - v 2 t)z) 

[|ri±| 2 + 71(2 - Uli) 2 ] 3/2 [l r 2 ± I 2 + 72 ( z - V 21 ) 2 ] 3/2 

B = p! x E x + (3 2 X E 2 (5) 

For ultra relativistic bunches, /3\ :2 —> 1, 77 _•> 1. The 

electric and magnetic fields take on the form of sharply peaked 
transverse pulses centered around the z position of the volume 
elements. The magnetic fields become of the same magnitude 
as, and perpendicular to, the electric fields. 

Physics processes occur as the bunches collide and volume 
elements V\ and V 2 overlap. The true fields that exist during 
the bunch collision are complicated, as the bunches electro¬ 
magnetic interactions produce field distortions. Analytic calcu¬ 
lations of transition rates that take into account these true fields 
exactly, are unfeasible. The calculation strategy is to determine 
transition rates of physics processes of interest within each vol¬ 
ume element. To then calculate the transition rate across the 
whole bunch collision, a simple sum is required for all the vol¬ 
ume elements encompassing the bunch collision. 

This strategy allows the electromagnetic fields within the 
volume elements to be approximated. If the volume elements 
are small and relativistic enough, the field’s longitudinal com¬ 
ponents can be neglected and the field strengths can be con¬ 
sidered constant. Such approximations allow the electromag¬ 
netic fields of each volume element to be described by a con¬ 
stant crossed field. The 4-potential representing two such over¬ 
lapping volume elements is a sum of two terms. Each term 
contains information about the propagation direction of each 
bunch, which can vary independently depending on the chosen 
reference frame 


A® = ea^h-x) + ea 2 p(k 2 -x) (6) 

aip = (0, ai), a 2fl = (0,a 2 ) 
kip = (wi,fci), k 2 p = (W2, k 2 ) 

The bunch collision may not necessarily be centred with re¬ 
spect to the centre of each bunches transverse cross-section. So, 
the three potential of the field of each volume element may have 
any direction as long as its transverse to the propagation direc¬ 
tion of its own bunch. The 4-vectors ki, k 2 are null 4-vectors 
whose 3-vector lies along each bunch’s propagation direction. 

In the calculation of Lorentz invariant transition rates of 
processes embedded in the fields, scalar products will be formed 
from the electromagnetic field tensor. These typically appear as 
a dimensionless parameter, 

T = sJ{eFp vVv Ylm z (7) 

For the constant crossed fields considered in this paper (equa¬ 
tion | 6 j, the invariant quantitiy T is, 

™ fc i ■ P + a 2p k 2 ■ p) 2 

i = e - 5 - (o) 

In the interest of precision, and since the electromagnetic 
fields of highly compressed charge bunches are intense, it is 
desirable to adopt a QFT which takes into account the charge 
bunch fields exactly. The most useful QFT to use is the Furry 
picture which is non-perturbative with respect to the external 4- 
potential. The Furry picture will be used to generate equations 
of motion for spin 1/2 particles. 

The solutions obtained from the equations of motion will 
then be applied to the particular case of two constant crossed 
electromagnetic fields associated with ultra-relativistic charge 
bunch collisions. The solutions will be used to calculate a new 
analytic expression for the beamstrahlung transition rate in both 
charge bunches, which reduces to the regular expression in ap¬ 
propriate limits. 


3 Dirac equation solutions in collinear fields 

In the Furry picture l4 HT7llT8l|T9l , the electromagnetic gauge 
field is separated into a part Ap that interacts perturbatively 
with coupled fields, and an external part A® whose interaction 
with coupled fields is calculated exactly. Thus, the Fagrangian 
density for a charged spinor of mass and charge m, e is lf 20 l . 

^dirac = $(*$ - eX - rn)ip (9) 

from which the Euler-Fagrange equation gives the equation 
of motion, 

(i$ — eJ/i — m) t/^pp = 0 (10) 

For an external potential described by a single plane wave, 
the solution of the minimally coupled Dirac equation (equation 
[Tol l is the well known Volkov solution lf2~il . <Zv, which differs 
from the usual free fermion solutions by the inclusion of an 
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extra phase and extra spinor E p . For an electron of momentum 
Pp — (e> P) i n an external field of 4-momentum k p = (w, k) and 
with uf p being the positive and negative energy Dirac spinors 
with spin r, the positive and negative energy Volkov solutions 
are. 


= / 


■Ep{4>) ~rp 


ut. e Tip ' x 


e£(4>) = 


it 


2 (k-p) 


2e(27r) 3 

4 , 2eA e (£)'P =f e 2 A e (£) 2 


T 1 


r 


2 k-p 


dZ 


(ii) 

= k-x 


The fermion solutions in two co-propagating plane wave 
fields, V'ocp- were obtained for the special case of orthogonal 3- 
potentials mm and applied to various combinations of laser 
fields. The solution obtained contains a product of Volkov E p 
functions. 


Ap — Aip^i) + Ai • A 2 — k\-k2 — 0 


V’oCP^l’ •fo) — 


)l,2 = ki^-x 

E %(<t>A E t v (<t>A u± p 




sJ2e{2-k) 3 


( 12 ) 


EtMi) = 


IT 


eAlf 

2 (ki-p) 


r4>i 2 eA;(£;) - pTe 2 ^;) 2 

^■1 - ^k—p - dii 


However, mm were overly constrictive, since the re¬ 
quirement of orthogonal 3-potentials is not necessary. The 4- 
potential of two (and indeed any number of) collinear plane- 
wave fields of 4-momenta kip, k 2 p can be rewritten as a func¬ 
tion of a null 4-vector n ; , whose 3-vector part is parallel to the 
common propagation direction of the collinear fields, 

kip = wi rip, k 2p = w 2 rip, rip = (l,n), re 2 = 1 

Mp{4>i) T A 2 p{<t> 2 ) = A e p (n-x) (13) 


One then posits that the general positive and negative en¬ 
ergy solutions ipQ P of the Dirac equation minimally coupled to 
the combination of two co-propagating external fields A e , are 
Volkov solutions with respect to the 4-vector rip. 


V'cp (n-x) = 


j± = 
^ px — 


it 


A e {n-x)if 


m T± ,,± pTip-x 
2e(2n) 3 P x r P 

*2 A'K).pTA'($f 


(14) 


n-p 


Ti/" 


2 n-p 


d£ 


An explicit proof that this is indeed the case for two collinear 


plane wave fields of general form, is presented in Appendix A 


To use this solution in the calculation of general particle 
physics processes in the fields of the two overlapping bunches, 
canonical commutation relations and quantisation are achieved 
within the Furry picture in the usual fashion li20l . A prerequi¬ 
site for the quantisation is the requirement of orthogonality and 
completeness of the solutions tpcp- 

Il22ll proved the orthogonality of the Volkov solutions for 
fermions in a single plane wave electromagnetic field on the 


equal time hyperplane. The proof proceeded by finding a new 
basis for the Volkov spinor which resulted in the Volkov phase 
when acted on by its hermitian conjugate. l23l achieved the 
same result by rewriting the standard Volkov spinor in terms of 
the derivative of the Volkov phase. Completeness of the Volkov 
solutions on equal time hyperplanes was proven by EH- 


Also required for the calculation of transition rates, is an 
explicit form of the general solution for the constant crossed 
fields expressed in equation[6] The explicit solutions are 


ij.’cp(n-x) = 


_ ™ _ 7± w ± e Ti(p-x) ns) 

2e(27r)3 P x r P V ’ 


J± = 
^ px - 




it 


(w\f 1 +o; 2 ^ 2 ) t/i (n-x) 


n-p 

(uJiai-p + CD 2 G 2 -p) (n-x) 2 q= + UJ 2 CL 2 ) 2 (n-x) 3 

2 n-p 


In the next section, the explicit fermion solutions ip^p will 
be applied to the beamstrahlung process, which is a 1-vertex 
processes within the Furry picture. The beamstrahlung transi¬ 
tion rate for two bunch fields will be compared with that for 
one bunch field. 


4 The beamstrahlung in overlapping charge 
bunch fields 

The beamstrahlung (figure [3]) is the photon radiation from in¬ 
dividual charged particles at a collider interaction point due to 
interaction of the charged particles with the fields from the col¬ 
liding bunches. Usually, the field of the particle’s own bunch 
is neglected and only the interaction with the oncoming bunch 
field is considered. Here, the beamstrahlung transition rate in 
both bunch fields will be calculated with the aid of the fermion 
solutions obtained in section [3] 



Fig. 3. The beamstrahlung Feynman diagram. 


The matrix element of the beamstrahlung process can be 
written down with the aid of the Feynman diagram (figure [3]), 

M fi = -ie f d 4 x $+ v {x,p f ) 1 pil } +,{ X ,p i ) (16) 
































A. Hartin: Furry picture transition rates in the intense fields at a lepton collider interaction point 


5 


One could proceed, as usual, in writing the square of the 
matrix element and integrating over space-time and wave packet 
momentum. However, much of the work has been done for this 
process, where the transition probability was obtained for any 
number of overlapping external fields with a common propaga¬ 
tion direction ll25l . Writing the transition propbability as a rate, 
adapting the result for the 4-vector n ; ,, and allowing for a rel¬ 
ativistic initial particle expressed by 7 = 6,/ to , the transition 
rate per unit time is, 


W^beam = ~ 


47T7 

2 - 


d(j> w 
2 ?tL 


du 


—ttL 
l + (l+u) 2 
2(1 + u) 


(1 + «) 2 / 


iird\ 

2A 


F\ (A 


.exp 1 


2 (n ■ Pi) 


A- 


Fb(A + (j > w , 4> w ) 2 


+ ^3(A + ^>iu, 



The functions F\. F 2 , F 3 depend on the external field 4- 
potential through the relations. 


Fifa,^) = —^ [A e (^ v ) - 


F‘2 {tpv 5 ) — 




(18) 


is the centre of mass frame of the colliding electron-positron 
bunches. 


The Xi, T> parameters, one for each of the charge bunches 
(see figure [4]», are effectively ratios of the field strength expe¬ 
rienced by an individual charge, to the Schwinger critical field 
strength of E CI = 1.32 x 10 18 V/m. 


With 7/7 .2 being the relativistically boosted field strengths 
of the bunches in the centre of mass frame, 71^2 being unit 
null 4-vectors representing the propagation direction of each 
bunch, and the unit 4-momentum of an individual charge pi = 
(1 


X.2 = 


e|ai|(fci ■ pi) _ 7 ^ 1,2 


(Xl,2 • Pi) 


( 21 ) 



F 3 (^,0 w ) = -4 

m J&-W 


Fig. 4. Beamstrahlung process in two overlapping charge bunch 
fields 


where the 4-potential itself is a sum of the constant crossed 
fields of two ultra-relativistic charge bunches and is a function 
of the unit vector scalar product n.x via the relation, 

= (wiai^ +u} 2 a 2 p) {n ■ x) (19) 


Typically for planned future linear colliders such as the In¬ 
ternational Linear Collider f26l . 7-E/ 2 /£/ r is of the order 0.1. 
So 'yEi^/E CI = 0.1 is used here in a numerical consideration 
of the beamstrahlung transition rate. 


Substitution of all relevant relations into the expression for 
the transition probability, and integration over the (f> w , A pa¬ 
rameters, results in the final expression. 


llbCLlTTl - 



r°° du 

Jo A + ^ 



1 + (l + tx) 2 d 
z(l + u) dz 


Ai (z) 


where 


(Tia! +T 2 a 2 ) 2 


1/3 


( 20 ) 


In ITbeam, the argument of the Airy function, z, is the only 
object that depends on the field strengths of the charge bunch 
fields. These field strengths are expressed by the parameters 
Ti, T‘j, where a l5 a 2 are unit 3-vectors lying along each bunch 
field’s 3-potential. 


The combination (Tiai + T 2 a 2 ) 2 is the Lorentz invariant 
quantity, T, of equation [8] IT'bcam reduces to the usual beam¬ 
strahlung expression when the field of the initial particle’s charge 
bunch is neglected (Xi —» 0). 


Since the transition rate is a Lorentz invariant quantity, the 
final analytic form of Wbeam is the same in all reference frames. 
The reference frame usually considered for collider processes. 


If allowance is made for an angular variation of an individ¬ 
ual radiating particle with respect to the propagation direction 
of its own bunches (9 in figure[4|, it is possible to produce a nu¬ 
merical comparison of the beamstrahlung due to the oncoming 
bunch only with the beamstrahlung that takes into account both 
charge bunch fields. Figure [5] shows the beamstrahlung transi¬ 
tion probability as the angle 9 varies. 

A full analysis of the significance of the new analytic form 
for the beamstrahlung transition rate is to be left to future work. 
Nonetheless the aim of this section, that of showing that fermion 
solutions in two overlapping charge bunch fields lead to ana¬ 
lytic and numeric differences in transition rates, has been achieved. 
Moreover, the use of fermion solutions in two charge bunch 
fields within the Furry picture, provides a methodology for cal¬ 
culating the transition rate of other collider processes in the 
presence of the intense electromagnetic fields of colliding charge 
bunches. 


5 Conclusions 

The field strengths of the charge bunches of future lepton col¬ 
liders are strong enough to potentially affect all collider pro- 
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Fig. 5. beamstrahlung transition rate per unit time due to initial 
particle angle for 250 GeV electrons and 7 ^ 1 . 2 /A’cr = 0.1 


cesses taking place in their presence. This has prompted a de¬ 
sire to theoretically treat particle physics processes that occur 
in the overlapping bunch fields at the interaction point, non- 
perturbatively with respect to the charge bunch fields. Such a 
treatment is provided by the Furry picture, which calculates the 
exact contribution from the fields. The Furry picture reduces to 
the normal interaction picture in the limiting case of vanishing 
external field strength. 

Employment of the Furry picture at collider interaction points, 
requires solutions of the field equations of charged particles 
coupled to both charge bunch fields. Since these fields are tied 
to massive sources, it is possible to choose a reference frame in 
which the charge bunches co-propagate and the problem sim¬ 
plifies. Solutions of the equations of motion for spin 1/2 par¬ 
ticles in two co-propagating fields with co-planar 3-potentials 
were obtained. In the limiting case of orthogonal 3-potentials 
and the vanishing of one of the fields, the general solution re¬ 
duces to previously obtained solutions illlll2ll2Tl . 

The fermion wavefunction in two co-propagating fields is 
written in terms of a null 4-vector whose 3-vector part has a 
unit norm and lies along the common propagation direction 
(equation p~4]). These solutions are orthogonal and complete on 
equal time hyperplanes. Canonical commutation relations and 
standard quantisation were assumed to follow from the orthog¬ 
onality and completeness properties. 

The solutions were used to calculate a more general form 
of the beamstrahlung transition rate for an electron radiating a 
photon in the fields of both colliding charge bunches. The an¬ 
alytic difference between the one charge bunch beamstrahlung 
and two charge bunch beamstrahlung is encapsulated by the 
Lorentz invariant parameter T, which is a function of the elec¬ 
tric field strengths of the charge bunch fields and the angle be¬ 
tween them. 

A numerical comparison of the beamstrahlung transition 
rate in one charge bunch field, to that in two charge bunch 
fields, showed differences for beam parameters expected at a 
future linear collider. This difference, depending on an angle 
between the propagation direction of the initial particle and that 
of its charge bunch, may be important for precision calculation 


required at linear colliders. A future detailed study will exam¬ 
ine this process and its ramifications in greater detail. 

A crucial question, given the precision physics program 
planned for future linear colliders, is, what is the effect of the 
strong charge bunch fields on collider process in general? Build¬ 
ing upon the work done here, it will be natural to consider tran¬ 
sition rates of general physics processes using solutions of the 
equations of motion in collinear fields within the Furry interac¬ 
tion picture. 

Apart from the beamstrahlung, the coherent pair produc¬ 
tion, photon absorption and decays in final states can all be con¬ 
sidered first order Furry picture processes. Second order pro¬ 
cesses involving an incoming electron and positron in the initial 
state, will necessarily interact with both charge bunch fields. 
The fermion solutions in two charge bunch fields used here, 
will be a prerequisite for the calculation of transition rates. 

For ongoing work, the wavefunctions of charged spin 0 and 
spin 1 particles in two charge bunch fields will have to be spec¬ 
ified. Two point Furry picture correlation functions leading to 
the particle propagators in two charge bunch fields are also re¬ 
quired, as are techniques needed to calculate higher order Furry 
picture transition rates. This program of work will be embarked 
upon in forthcoming papers. 


This work was supported in part by DFG Sonderforschungsbereich 
SFB 676 


Appendix A Solution of the Dirac equation 
in two collinear fields 


Required, is the wavefunction solution, t/>cp, for a fermion embedded 
in two co-propagating, plane-wave, electromagnetic fields of general 
form. Linear combinations of the 4-potentials can be formed in order 
to achieve orthogonal components. With general functions fi, / 2 


A^ = 01^/1(01) + a 2M / 2 (> 2 ) = Ai(</>i) + A 2 (</> 2) 

= aiM^il/iC^i) + a2 M |a 2 |/ 2 (</> 2 ) 

= a +ll f+{(/)!, cf) 2 ) + a-,,cf)2) ( 22 ) 


where a± M = 


ai^ ± a 2 ^ 

V2 


f± = 


a±-fci,2 = a+-a_ = fci-fe = 0, 

CL ^ — CL 2 — 1 


I q i l/i ± \a2\f2 

s/2 

01,2 = kl2'X 


If ipcp is a solution of the Dirac equation minimally coupled to the 
external field, it is also a solution of the second order Dirac equation, 

[(*# - eA e ) 2 + m 2 + ip C p = 0 

with = ^( 7 ^ 7 ,, - 7 „ 7 P ) 

The reformed external field 4-potential and the orthogonal rela¬ 
tions of equation[ 22 ]allow separable equations to be formed, 

[(*# - ea±f±) 2 + m 2 + ie<jbf±<ji ± j t/> C p = 0 (23) 
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There are two sets of solutions corresponding to positive and neg¬ 
ative energy solutions. Proposing i/j„ as a product of the free, fermion, 
plane wave solutions and general functions of the external field, G± {(j> i ; 

i } cp{ < t , u fa) = e Ti [ p ' x+G+(,il ’ 02)+G - Wl ’ 02) ] (24) 

substitution into the second order Dirac equation yields first order 
partial differential equations, 

[2 (fa -p) JL + 2 (fa -p) ^-] G± (0! ,fa) = F ± (fc , </> 2 ) 


where F± = 2e ( a±-p ) f± — e 2 a±f± + ie$f±tf 


(25) 


Choosing a structure for the unknown functions G ± that matches 
that of the known functions F±, separate differential equations can be 
formed, 

G±(<j>i, fa) = g±i(<j>i) T g±2 (02) + g± 3 (fa, §2) 

dg± 1 2 e(a±-p) |ai|/i - e 2 a^.\a 1 \ 2 f 2 + ief 1 \a 1 \ft 1 ^, ± 


d<t> 1 2 v / 2 (fci-p) 

dg ±2 ±2e (q± -p) |a 2 |/ 2 - e 2 a^_ |a 2 1 2 /| ± ie/ 2 |a 2 |#6 2 #4 ± 


d<t> 2 


2 % / 2 (fe 2 -p) 


(26) 


r o ,, \ d ,0/7 X d 1 e 2 a||a 1 ||o 2 | 

[ 2 ^ ^ + 2 (fc2 ' P) 5±3 = T ---* /2 


/2 


These differential equations are easily solved. Forming the sum 
G+ T G-, the positive and negative Dirac wavefunction solutions in 
two co-propagating plane wave electromagnetic fields of general form 
are. 


V’cpW’Ij fa) = cEt p (fa)Et p (fa) ut p e^- e2 ^- a2H ] (27) 
IT 


Ef p (fa) = 


2(fej-p) 




2ki ■ p 


Hs j + 

Limiting cases of the above general solution were obtained earlier 
in the literature. Restricting the 3-potentials to orthogonality ai-a 2 = 
0, the term containing the H function drops out and a product of 
Volkov solutions is obtained mm. Additionally, a single Volkov 
solution results from removing one of the fields by setting |ai | or 1 0,2 1 
to zero. 

For collinear fields, the energy can be extracted from photon 4- 
momenta leaving a null 4-vector n M = (1, n), whose 3-part is parallel 
to the collinear direction 


fa ,2 p=wi ,2 <£ 1,2 =wi ,2 n-x, fci, 2 -p=wi ,2 n-p 


(28) 


In equation[27] making shifts of integration variables, —>■ , \ 

ojix and gathering terms, the general solutions reduce to a Volkov so¬ 
lution in terms of the 4-vector n u . 
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